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ABSTRACT 

We describe a combined dynamic atmosphere and maser propagation model of SiO 
maser emission in Mir a variables. This model rectifies many of the defects of an ear- 
lier model of this type, particularly in relation to the infra-red (IR) radiation field 
generated by dust and various wavelength-dependent, optically thick layers. Modelled 
masers form in rings with radii consistent with those found in VLBI observations and 
with earlier models. This agreement requires the adoption of a radio photosphere of 
radius approximately twice that of the stellar photosphere, in agreement with obser- 
vations. A radio photosphere of this size renders invisible certain maser sites with 
high amplification at low radii, and conceals high-velocity shocks, which are absent in 
radio continuum observations. The SiO masers are brightest at an optical phase of 0.1 
to 0.25, which is consistent with observed phase-lags. Dust can have both mild and 
profound effects on the maser emission. Maser rings, a shock and the optically thick 
layer in the SiO pumping band at S.lSiitmappear to be closely associated in three out 
of four phase samples. 
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1 INTRODUCTION 

We present preliminary results from a significantly im- 
proved model of SiO maser emission from the circumstel- 
lar envelopes of Mira-type, and similar, long-period vari- 
able stars. The predecessor of this model, hereafter the 'old 
model', was used to compute synthetic spectra and images 
of SiO maser emission from a model star with parameters 
based on those of o Get ("Grav etal ] ll995al :[ Humphrev s et al.l 
[l996l). With additional phase information, the variation of 
the spectra a nd maps was computed over the stellar pe- 
riod of 332 d (iHumphrevs et al.lll997al : iGrav Humphrevsl 
|200Q| : iHumphrevs et al.ll2002l ). The hydrodynamic solutions 
which provided the input physical conditions for the maser 
sites in t he old model were taken from envelope computa- 
tions b^Bowen ( 1988r), except in the earliest of these pa- 
pers, I Grav et al.l ( 1995al l , which used similar calculations by 
IWillsonI (|l987f r 

The old model used the (radial) density, velocity, and 
gas kine tic temperatu re p rofiles from hy drodynamic solu- 
tions bv lBowenI (|l988l ) and lWillsonI (|l987l ). These hydrody- 
namic solutions also provided a radiative equilibrium tem- 



perature. Teg, on the basis of geometrical dilution of the 
stellar continuum, and on the opacity of the gas, but no es- 
timate of the abundance of SiO. The radius of the stellar 
photosphere was recalculated at each time step, and defined 
as the radius where the optical depth equalled 2/3, but the 
eflFective temperature of the photospheric layer, whatever 
its radius, had a constant value, independent of the stel- 
lar phase. Dust was included in the pulsation model, and 
was dyna mically s ignificant at radii beyond the condensa- 
tion zone (|Boweni ri988). At radii where dust was present, 
the dust temperature was set equal to Teq. The eflFect of dust 
was exerted through a spectrally averaged radiation pressure 
cross section, a function of Teq and the dust condensation 
temperature. The solutions did not specify a wavelength de- 
pendence for the dust opacity, or a dust number density, 
though this latter quantity could have been calculated, us- 
ing sensible assumptions for the dust mass fraction, grain 
radius and mineral density (see Bow en|[l988l V Later Bowen 
models impro ved upon the version discuss ed above in sev- 
eral respects (jBowenl [T990l : IWillson II2OOOI ). but these later 
models were never used by us for maser calculations. 
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The old model made several simplifications and approx- 
imations, not inherent in the Willson and Bowen solutions, 
for constructing the IR radiation field for maser pumping. 
The stellar component of IR radiation field assumed a fixed- 
radius stellar photosphere, and was modelled as a black- 
body at the (constant) stellar temperature of 3002 K. Opti- 
cally thin dust emission was crudely modelled as a black- 
body at a single temperature, weighted by a power-law 
in wavelength (Humphrey s et a l. 1996). As a consequence 
of these approximations, the radiative pump rate in any 
SiO rovibrational transition had no intrinsic variation with 
phase. Although the collisional part of the maser pump was 
modelled more accurately than the radiative part, a radia- 
tive com ponent to the maser pumping was, however, always 
present (|HumDhrevs et al.lll99d l2002l ). The old model used 
a constant value of 10~^ for the abundance of SiO, inde- 
pendent of both radius and phase, and, for the purposes 
of pumping calculations, hydrogen was assumed to be in 
molecular form, as no suitable set of rate coefficients was 
then available for the H+SiO system. 

The stellar phase in the Bowen hydrodynamic solutions 
was defined a s zero at the time a new shock left th e stellar 
photosphere (jBowenlHissI : iHumphrevs et al.lll996l \ and it 
was never possible reliably to relate this model phase to 
the optical, or near IR, light curves of the star, despit e 
several attempts t o do so, for example iGrav et al I (|l998h . 
IHumphrevs et al" It was therefore impossible to an- 

swer questions related to observed links between the light- 
curves of Si O masers and the optical and near-IR continuum, 
for example IPardo et al.l (2004). Polarization and the pres- 
ence of minor isotopomers of SiO were also not included in 
the old model; these aspects of the problem remain absent 
in the new model, and will not be considered further in this 
work. 

Despite the drawbacks listed above, the old model 
was surprisingly successful when co mpared with obser- 
vations. The single phase version (jGrav et all Il995al : 
IHumphrevs et al] 1 19961 ) predicted new masers with up- 
per rotational states above J=6. Such masers were sub- 
sequently detected in a number of stars, in cluding ex- 
ample s from Miras and semi-repular variables ( Grav et aP 
1 1995bl : IHumphrevs et al.lll997bl : Icrav et al.lll999l l. The ring 
structure of the SiO masers , found frequently in VLB I 
observations ([Diamond et al.l 1 19941 : iGreenhill e t al.l I1995I : 
iBoboltz et al.lll997l ). which indicates predominantly tangen- 
tial amplification from a restricted range of radii, was also 
strongly supportive of the model. It is worth noting here 
that the model generated potential maser sites randomly 
throughout a sphere of significantly larger radius than the 
ring in which masers appear. If the model had been seri- 
ously inappropriate, it could have produced, for example, 
a smooth radial distribution, or, for predominantly radial 
amplification, a concentration of bright maser spots directly 
over the central star. The ring radius for the v — 2^ J — 1 — 
masers was found to be smaller, in a number of stars, 
than for v — 1 , J = 1 — , in agre ement with the model 
dPesmurs et al." 1999, 2000; Yi et al.l l20Q0l : ISoria-Ruiz et al.1 
12004 : Boboltz k Wittkowski 2005). The model prediction 
for the ring radius of the masers in the v — l^J — 2 — 1 line 
at 86 GHz, is that it should have a comparable, or slightly 
larger radius than the ring formed by the v — l^J— 1 — 
masers, and be significantly less well populated. The sec- 



ond prediction seems t o be in ageement with observations 
(ISoria-Ruiz et al.1 l2004l) but not the first. Observations of 
R Leo (iSoria-Ruiz et al ] 12007*) confirm that, in this star at 
least, the v — l^J — 2 — 1 transition forms in a ring further 
from the star than the J = 1 — transition from the same vi- 
brational state. Single-dish observati ons of a number of Mi- 
ras w ith different dust temperatures (|Nakashima k, Deguchil 
I2OO7I ) showed a correlation of the flux in the 8 /xm continuum 
and the J = 1 — masers from v — 1 and 2, which strongly 
suggests a radiative component to the maser pump. 

The maser line ratios appear to be acceptable, given 
the uncertainties in what constitutes the background for 
unsaturated masers, except for the v = 2,J = 2 — 1 line 
which is much too strong when compared against observa- 
tions. However, the weakness of this maser in real sources is 
very probably due to a line overlap with a water transition, 
for example Soria- R uiz et al.1 ( |2004l ). and a molecular model 
of water was (and is) not included in our maser model. 

When phase dependence is added, the old model re- 
mains broadly in agreement with interferometric observa- 
tions, though with some discrepancies. For example the 
model prediction of smaller maser rings, att' = 2,J=l — 0, 
compared with ^; = 1, J = 1 — 0, over the stellar cycle, is in 
agreement with observations (ICotton et aD 20041 : lYi et al.l 
I2OO5I : ICotton et al. I I2OO6I : IWittkowski et al.ll2007l ). as is the 
prediction of an r.m.s. amplitude variation of 5-10 per cent 
in the maser ring radius. These maser rings are also found to 
be inside the inner radius of the silicate du st shell, but AI2O3 
dust may be co-located with the masers (|Wittkowski et al] 
2007). However, the predicted relation between SiO maser 
ring radius and SiO maser luminosity is not well reproduced, 
and varies both between differe nt stars, and for differen t 
pulsation cycles of the same star (|Cotton et alJl2003j2006l V 
Long term monitoring of the v — 1, J — 1 — line in the Mira 
TX Cam (Diamond & Kemball liooi ) shows motions domi- 
nated by shock-driven outflow for part of the stellar period, 
whilst the remainder is dominated by infall under gravity. 
This behaviour agrees with the model assuming a reason- 
able link between the model and optical phases, though the 
variation in the maser ring radius is somewhat larger (18 per 
cent) in TX Cam. It should also be noted that TX Cam ex- 
hibits considerable asymmetric behaviour, and some maser 
features move outwards, for example, even during the infall 
do minated phase. If the model to optical phase relation used 
bv lYi et "all (|2005l l is adopted (model phase = 0.78 at opti- 
cal maximum) then the disrupted, weak maser ring structure 
found near optical minimum in TX Cam is also consistent 
with observations, though any linking of the old model and 
optical phases should be treated with great caution for rea- 
sons discussed above. 

Single dish monitoring observations suggest that shock 
waves, at least of a type which take a significant fraction 
of the period to cross the gap between the stellar photo- 
sphere and the maser zone, cannot be responsible for the 
maser pumping, though some observations still suggest a 
coupUng bet ween the pump and t he kinematics of the enve- 
lope gas. See Cotton et all (|2006l l for a detailed discussion 
of the problems involved in relating results from interfer- 
ometric and single-dish observations. Whilst it is not true 
to say that the old model is entirely collisionally pumped, 
a major drawback is that it cannot explain correlation be- 
tween the light curves of SiO masers and radiation in the 
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IR pumping bands . It is clear that such correlations exist: 
I Par do et alj (|2004h found that the optical, near-IR contin- 
uum and SiO masers share the same period in a sample of 
13 Mira variables. They also found the IR and SiO maser 
peak luminositi es coincided t o with in 0.05 of a stellar period. 
The IR data in lPardo etaP (|2004l l is at 4.9 and 3.5 /xm, and 
therefore does n ot include the cru cial = 1 pumping band 
of SiO at 8 /im. 'Mcl ntoshl (|2006al l finds that the maxima in 
the 43 GHz masers from v — 0, 1, 2 occur simultaneously to 
within 0.066 periods in o Get. As the v — maser ring was 
apparently much larger than those dX v — 1 and 2, this sug- 
gests that the pumping is predominantly radiative. However, 
the old model makes no prediction about v — masers, and 
the difference between the ring radii v — 1 and 2 is suf- 
ficiently small (^0.1 to ^0.2 AU, dependent on phase) that 
a shock of reasonable speed could propagate between them 
in 0.066 of a period (22 d). These observations (Mcintosh 
|2006a|) also did not detect any disruption of the maser ac- 
tivity, due to shock arrival, near optical minimum. jFurther 
analysis with the same dataset from o Get (|McIntoshll2006bl ) 
finds that SiO maser flux density in the v — 1^ J — l — tran- 
sition is co rrelated with the velocity centroid of the emission. 
iMcIntoshI (|2006b :> offers two possibilities for this correlation: 
either the pump affects the kinematics of the maser region 
(that is, shocks strongly influence maser amplification, prin- 
cipally through disruption of velocity-coherent gain paths, 
but also by heating and compressing the gas) or that the 
maser pump and (independent?) kinematic effects combine 
to produce the observed maser spectra. The current authors 
find this second possibility hard to reconcile with the model, 
since in the large velocity gradient (LVG) approximation 
used, the Sobolev optical depth in any direction is depen- 
dent upon the velocity gradient, and therefore the kinemat- 
ics are intimately coupled to the pumping scheme. Similar 
observations for the Mira R Gas ([Mcintosh k Havnesll2007l ) 
largely agree with those from o Get. However, it is interest- 
ing that for one maximum, the peaks of the v — 1 and 2, 
J = 1 — masers are significantly separated in phase. 

The pulsating atmospheres of Miras are very ex- 
tended because of dynamic events, such as shocks. These 
extended atmospheres exhibit phase-dependent tempera- 
ture and density stratifications, optically thick roolecu- 
lar layers, and dust formation; see for example ISchold 
(^03). Observed radii of Miras have been found to 
depend on both phase and wavelength of observation, 
which is consistent with the predictions of dynamic mod- 
els, in which molecular layers lie abo ve the c ontinuum- 
forming laye r ([T hompson et al. 2002: Perrin et al. 2004; 
IOhnakal[2003 : lohnaka et al.ll2(306bl: IWittkowski et al.ll2008l) 
In a series of experiments (iBoboltz &: Wittkowski| l2005l : 



a series ^ 

IWittkowski et al.1 boOTi ) mid-IR interferometric measure- 
ments made with the VLTI were made near contempora- 
neously with VLBA observations of the 43 GHz masers in 
V — 1 and 2 towards the Mira S Ori. Models (jlr eland et al.l 
l2004al lbh were fitted to four observing epochs, with an esti- 
mated error of <0.15 pulsation periods between the model 
and optical phases. Both the IR photosphere radius, defined 
as the radius at 1.04 /xm, and the radius of the optically 
thick layer near 8 /xm (the A^; = 1 SiO pumping band) vary 
in phase with the optical light curve, as far as can be judged 
from the limited phase sampling. Importantly, the radius 
of the 43 GHz maser shells is slightly larger than this 8 /xm 



quasi-photospheric band radius at all four observing epochs. 
If this optically thick layer represents high-density post- 
shock gas, then the position of the shock would be tightly 
correlated with the light curve at this frequency, very possi- 
bly leading to a resolution o f the radiative /c o llision al pump- 
ing debate. Recent work bv^ Woodruff et aP (|2008l l suggests 
that, with certain caveats, the conclusions obtained for S Ori 
can be extended to a sample of other Mira variables. 



2 IMPROVEMENTS TO THE MODEL 

Our improved model removes many of the defects inher- 
ent in the old version because of its restricted information 
about conditions in the circumstellar envelope. In particu- 
lar, the new version includes radii for several optically thick 
layers, for both continua and restricted- wavelength bands, 
as a function of phase. The set of bands includes all the IR 
pumping bands of SiO which are accessible to the model. 
Other features of the new model include an accurate esti- 
mate of the IR radiation field from dust, sufficient chem- 
istry to provide abundances for SiO and its main collision 
partners, and new rate coefficients for collisions of SiO with 
hydrogen atoms. A link between the model phase (in terms 
of the shock position) and the stellar light curve is implicitly 
included in the model data. 

At this point we digress to make some definitions re- 
garding the optically thick layers introduced above. We re- 
serve the word 'photosphere' to mean a continuum-forming 
layer, and we always use this word with a modifier which 
specifies the region of the spectrum to which it applies. In 
particular, we define 'IR photosphere' to be the continuum- 
forming layer at 1.04 |j,m where the optical depth is equal to 
1. This wavelength was chosen because it is in a region little 
contaminated by molecular line absorption. We use 'radio 
photosphere' to be the continuum- forming layer in the radio 
region. This is less well defined, but, in the context of SiO, 
can be taken to mean the region of the electromagnetic spec- 
trum occupied by the pure rotational lines of SiO in general, 
and the maser lines at 43 and 86 GHz in particular. When 
referring to other work, for example the hydrodynamic so- 
lutions of the old model, which does not match the precise 
definition of the IR photosphere, we use the expression, 'stel- 
lar photosphere', to describe the continuum- forming layer in 
the optical and near IR. Where photospheric radii have al- 
ready been mentioned in the Introduction, they conform to 
the above definitions. 

Each SiO pumping band has a quasi-photospheric, op- 
tically thick, layer, due to the presence of large numbers of 
molecular lines, mostly from water. These layers have radii 
computed as r a = 1 radii based on the definition of the fil- 
ter radius in IScholz Takedal ((l987|) 

using a box filter of 
width 0.1 ^Lm. When referring to these optically thick lay- 
ers, the form of words we adopt is 'optically thick band ra- 
dius', or 'optically thick layer'. When referring only to the 
radii of such layers, we write a shorter form, specifying the 
wavelength of the band in microns, followed by 'radius'. For 
example, to refer to the radius of the optically thick molec- 
ular layer in the /\v — 1 band of SiO, we use the shorthand 
'8.14 iJ,m radius'. 
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Table 1. Parameters for the M-Series models (jlreland et al.l 
l2004bh . The radius Rp refers to the parent star (see text). The 
effective temperature and luminosity are computed for the non- 
pulsating parent star. 



Parameter 


Symbol 


Value 


Period 


P 


332 d 


Mass 


M/Mq 


1.2 


Metallicity 


Z/Z0 


1.0 


Luminosity 


L/L0 


3470 


Radius 


Rp/R0 


260 


Temperature 




2750 K 



2.1 Hydrodynamic Solutions 

The hydrodynamic solutions u sed in the pre s ent wor k are 
based on those pub lished in llreland et al ] (|2004al ) and 
Ireland et al. 1 (|2004bl ). We are particularly interested in the 
M-series of models, drawn from the latter work. All these 
models have a pulsation period of 332 d, like the real star 
o Get and the old model. However, other parameters of the 
M-series models might be more appropriate for a model of 
R Leo. The abundance of SiO is computed by an additional 
code by Ohnaka. This code employs equilibrium chemistry, 
including 34 elements and 235 molecular species, and as- 
sumes solar abundances of the react ants. Optically thick 
band radii are provided in the SiO rovibrational pumping 
bands at 2.03, 2.74, 4.09 and 8.14 /xm, corresponding to vi- 
brational changes, respectively, of A^; = 4, 3, 2 and 1. 

The M- s eries o f models is decribed in detail in 
llreland et al.l (|2004b[ ) . The important stellar parameters are 
briefly listed here in Table [T] The radius and effective tem- 
perature refer to a non-pulsating 'parent' star, where these 
values are independent of phase. The parent parameters are 
based on a layer where the Rosseland mean optical depth 
is 1. There are 20 M-series models altogether, some with 
the subscript n, den oting a new model, first published in 
llreland et al. (2004b). Other mod els appear ed previously 
in either Hofmann et al. (1998) or iTei et al.l (|2003 ). Each 
model has a cycle (pulsation) number and a phase within 
that cycle. At present the samples range from phase 0.49 
of cycle to phase 0.50 of cycle 2, giving ten models per 
cycle, though these are not evenly spaced i n phas e . Mode l 
parameters are listed in Table 1 of Irel and et al.l (|2004bl ). 
Phase values are given in terms of the optical light-curve 
and are estimated accurate to about 0.1 pulsation periods 
absolutely, but more accurately relatively (0.01-0.02 peri- 
ods) . An estimate of the phase- variable IR photospheric ra- 
dius is given as Ri.qa for each model. The basic models are 
dust free, and extend radially to 5 Rp, where Rp is the radius 
of the parent star. Th i s is a considerably smaller radial cov- 
erage than in iBowenI (| 19881 ) , but this has no consequences 
if the water and OH maser zones are ignored. We also note 
that the models we use are developed from grey dynamic 
models which extend to considerably larger radii than 5Rp, 
and which establish the initial pressure stratification for the 
refined non-grey models. 

In this work, we use the same model s which were used 
in the fit to S Ori (|Wittkowski et al.ll2007i ). These are M21n, 
M22, M23n and M24n, which have respective optical phases 
of 0.10, 0.25, 0.30 and 0.4. The stellar parameters of S Ori 
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Figure 1. Radial velocity of model M21n as a function of radius, 
showing the presence of outward-moving shocks. The shape of the 
outer shock appears smoothed in the model for computational 
reasons (dashed line: without smoothing applied). 



differ from those used for the M series of models in pulsation 
period, radius and mass. The M-models were used as the 
best available option to describe the dynamic atmosphere 
of S Ori, scaling the period to phases between and 1, and 
the radius to match the ob served angular diameter of S Ori. 
See also the discussion in IWittkowski et al ] (|2007h . We do 
not expect these scaling operations to lead to strong effects, 
but it should be kept in mind that remaining differences 
between model predictions and observations of S Ori might 
arise from this moderate mis-match of stellar parameters. 

As examples of the behaviour of the physical condi- 
tions in the model atmosphere we plot, for M21n, the ra- 
dial velocity (positive means outflow), kinetic temperature, 
and the number densities of SiO and its main collision part- 
ners as functions of radius in Figures [1] [2] and (3] respec- 
tively. A strong inner, and a weaker outer, shock are vis- 
ible in Figure Eat radii of 1.57 and 3.00 Rp (2.84x10^^ 
and 5.43x 10^^ cm). T he shock fronts are slightly smeared 
out in the models f see llreland et al. I l2004bl l and affect both 
the temperature and density graphs. The shape of the 
SiO number density curve follows closely that of H2. The 
optically-thick band radii in the SiO pumping bands at 2.03, 
2.71, 4.06 and 8.13 /xm are, respectively, (2.51,4.18,2.13 and 
3.98) X 10^^ cm, so the 8.13 /xm layer, at this phase, actually 
lies outside the radius of the strong, inner shock. 



2.2 Dust and the IR Radiation Field 

In I Wittkowski et"aL the dust-free hydrodynamic so- 

lutions described above were complemented by ad-hoc ra- 
diative transfer models of the d ust shell, using the ra dia- 
tive transfer code mcsim_mpi by lOhnakaet alT (|2006al l . A 
dust-free model of the M series was combined with a dust 
shell model using dust shell parameters so that the over- 
all models fitted the 4 epochs of observations of S Ori. 
The best-fitting overall models for the four epochs of ob- 
servation, A, B, C, D have parameters listed in Table O 
We therefore adopt as standard, the pairings epoch A+M22, 
epochB+M24n, epochC+M23n and epochD+M21n, and use 
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Figure 2. Kinetic temperature of the gas in model M21n: the 
constant temperature of 741 K at large radii is artificially set, 
since the equation of state is poorly known at the very low tem- 
peratures of the uppermost layers. 
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Figure 3. Number densities of H2 (solid line), He (dashed line), 
H atoms (chained line) and SiO (complex line) 



these in the work that follows, except in Section 1431 and Sec- 
tion [52] 

In the modelling of S Ori, combinations of AI2O3 and 
silicate dust were tested, but better fits were obtained from 
dust-shell models without silicates. Therefore, the epochA- 
D dust shells use AI2O3 dust only, which is assumed to 
be in the form of spherical grains of a single radius, a = 
0.1 /im. Each dust epoch provides, as a function of radius, 
the number density and temperature of the dust. Optical 
efficiencies for ab sorption and scattering were taken from 
I Ko ike et al." (1995") for wavelengths shortward of 7.8 /xm, and 
from Begemann et al. (1997) for longer wavelengths. These 
optical efficiencies, the table of number densities and tem- 
peratures, together with stellar parameters provide all the 
information necessary to compute the angle- averaged inten- 
sity of the IR dust radiation field in the circumstellar enve- 
lope (see Section |4|). 

Silicate dust can optionally be included in the dust shell 
calculations. When it is so include d, the optical effi c iencie s 
are the 'warm silicate' variety from lOssenkopf et al.l (|l992[ ). 
The silicate grains are assumed to be spheres with the 



Table 2. Standard pairings of dust epochs and hydrodynamic so- 
lutions with parameters of the model star and dust she lls. Obser- 
vational phases ((/)obs5 Column 2) refer to S Ori ( Wittko wski et aP 
^007). Theoretical phases (</>, Column 4) are the phases of the 
hydroynamic model, on the optical scale, fitted to each epoch. 
Column 5 is the ratio of the inner radius of the dust shell to the 
IR photospheric radius, and Column 6 is the visual optical depth 
of the dust shell. 



Epoch 


(pobs 


Model 




Rin/ Rl.04: 


TV 


A 


0.42 


M22 


0.25 


1.8 


2.5 


B 


0.55 


M24n 


0.40 


2.0 


2.5 


C 


1.16 


M23n 


0.30 


2.2 


1.5 


D 


1.27 


M21n 


0.10 


2.4 


1.5 



same single radius as their AI2O3 counterparts. Dust of this 
type has been used in a model of the envelope of the Mira 
GX Mon (Boboltz et al., in preparation). We note that the 
dust shell extends to substantially larger radii {Rin '^2Rp 
and Rout/Rin=1000) and, when present, the silicate dust is 
found at larger radii (inner radius '^4Rp) than the AI2O3 
dust. 

We plot, in Fig.|4l the number density of AI2O3 dust as 
a function of radii for the four epochs A-D. Fig. [5] shows the 
dust temperature as a function of radius for the same four 
models. 

As noted in the introduction, dust was dynamically im- 
portant in the hydrodynamic solutions of the old model, but 
not so with the dust shells attached to the dust-free hydro- 
dynamic solutions of th e current model. I n this connexion, 
we note that studies by ^ Ireland Scholj (|2006l l found that 
dust did not noticeably affect the density and temperature 
stratification of the model for radii < 2 Rp in Miras with 
parameters similar to o Get or R Leo. Prelimin ary results 
from very recent calculations ("Ireland et al.ll200 "8*) , which do 
include the formation and dynamical effects of dust, sug- 
gest that this view is correct, though only a small subset of 
the full phase and cycle data is currently available. These 
conclusions from models are suppo rted by measurements o f 
the inner dust radius with MIDI (^Wittkowski et al.| [2007l ). 
which place it in range 1.8-2.2 IR photospheric radii. 



2.3 Rate Coefficients for H+SiO 

For details of the ra te coefficients used for c ollisions between 
H2 and SiO, see Humphre vs et al.1 (|l996 h and references 
therein. These rate- coefficients were originally computed for 
He+SiO collisions, and can easily be reset to the latter sys- 
tem by changing the total and reduced masses. 

The inner regions of the circumstellar envelope contain 
substantial amounts of hydrogen in atomic form (see Fig.[3|). 
We have therefore included a set of rate coefficients for 
H+SiO collisions. For details of the computati ons of the po- 
tential surf ace and collis i on cro ss sections, see ljimeno et al.1 
(1999) and IPalov et all (|2002l ): for calcul ation of the rat e 
coefficients from these cross-sections see Palov et al. 
Rate coefficients can be computed for any rovibrational tran- 
sition in the first six vibrational states {v = — 5) for rota- 
tional states from J = — 40, and with a maximum change 
in rotational quantum number of |AJ| = 40. The accuracy 
of the rate coefficients degrades progressively above 5500 K, 
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Figure 4. Number density of AI2O3 dust as a function of radius 
for the four epochs A (sohd hue), B (dashed hue), C (chained 
hne), and D (complex hne) 
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Figure 5. As for Fig. ID but showing the dust temperature. 
Epoch B has the smallest inner radius, followed by Epoch A. 
Epoch A has the highest temperature at its inner radius. 

but this is not a problem with the hydrodynamic solutions 
used in the present work, where the maximum kinetic tem- 
perture is -6000 K (see Fig. [2]). The H+SiO rate coefficients 
are different, and systematically larger than for the He+SiO 
system. They are of similar magnitude for = tran- 
sitions, but larger by a factor of 10 to 100 for the vibra- 
tionally inelastic case. The expected decay in the magnitude 
of rate coefficients as Af increases is also rather slow, so that 
H+SiO collisions can efficiently make rather large changes 
in vibrational state. 



3 REMAINING PROBLEMS 

The improvements discussed in Section [2] do not solve all the 
problems with the model, and actually introduce one or two 
new, though minor, difficulties. Most of these relate to our 
lack of knowledege of the position of the radio photosphere. 
In addition, we still rely on the LVG approximation for the 
maser pumping code, and we have no rate coefficients which 
are based on the true H2+SiO collisional system, though 



some new calculations may soon become available ()Bienieckl 
I2OO8I ). 

3.1 The radio photosphere 

The new hydrodynamic solutions (Section 12. 1|) provide 
phase- dependent radii for the optically thick molecular lay- 
ers in the four rovibrational pumping bands which are pos- 
sible in our model, which has 200 levels distributed equally 
between the first five vibrational states (v — — 4). How- 
ever, the pure rotational transitions cannot be assigned to 
any of these bands, and the position of the optically thick 
layer, or radio photosphere, at the longer wavelengths (sub- 
mm to cm) of these rotational lines is much less certain. The 
radiation field from the dust can be ignored for the pure ro- 
tational lines because of the very low optical efficiencies at 
these wavelengths, but see Section [3.21 below. Therefore, we 
need to consider only the stellar radiation from the point of 
view of pure rotational radiative maser pumping. 

Estimates of the radius of the ra dio photosphere have 
been made by IReid MentenI (|l997l l , using three frequen- 
cies, all lower than any SiO pure rotational lines. The opac- 
ity for the radio photosphere is maintained by free-free in- 
teractions of electrons with neutral H and H2, the elec- 
trons originatin g from ionization of K and Na. Figure 10 of 
iReid &: MentenI (|l997f ) strongly suggests that the radius of 
the radio photosphere at the SiO frequencies (in the approx- 
imate range 40-1500 GHz) is smaller than at the measured 
frequencies, and may change by a factor of '^1-3 ov er this 
range of frequencies. Overall, iReid Mente^ (|l997l ) give a 
figure of 2 R* for the radius of the radio photosphere, where 
R* is the radius of the stellar photosphere. The figure c)f 2 R* 
still appears in more recent work ( Reid Menten 2007i V but 
the radio photosphere is also linked to the optically thick 
molecular layer which becomes opaque in the mid infra- 
red. With this limited information, it is difficult to make 
a good choice for the radius of the radio photosphere, but 
we have chosen to associate it with the molecular layer, as 
in iReid MentenI (|2007f ) , and we have set the radius of the 
radio photosphere equal to the optically thick band radius 
at the longest IR waveband (the 8 /xm radius) . This means 
that the radio photosphere, like the optically thick layers 
in the IR bands, moves with the stellar pulsations. We also 
note that if we can set ^ Rp, then this choice of radius 
in the radio region is in good agreement with the figure of 
2R*. However, given the uncertainties in the position of the 
radio photosphere, and the fact that the maser fines cover 
an order of magnitude in wavelength, so there is unlikely to 
be one optically thick radius for all these lines, we allow the 
choice of maser sites to vary from the smallest radius in the 
hydrodynamic solution, at each phase, up to the constant 
value of 5Rp. We can then investigate the effect of making 
the radius of the radio photosphere very small, and subject 
the site distribution to later exclusions on the basis of the 
larger 8 /xm radius discussed above. 

3.2 Maser backgrounds 

In the old model, the background amplified by the masers 
was set to a Planck function at the gas kinetic tempera- 
ture local to each maser site, without much justification for 
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this choice. However, given that the allowed masers in the 
model are a subset of the pure rotational transitions dis- 
cussed above in Section ITTI fthose with upper J ^ 10), the 
presence of a radio photosphere at radii similar to those of 
the maser zone makes this choice rather more acceptable, 
and we keep it for the new model. The alternative would 
be to integrate along the line-of-sight through the dust shell 
over all the dust and SiO spontaneous emission lying behind 
the maser site. However, this background is optically very 
thin, and would produce a far weaker background than opti- 
cally thick emission from the radio photosphere. Of course, 
for saturating masers, it does not matter what the back- 
ground radiation is. However, the amplifying column re- 
quired to reach saturation is significantly longer for a weak 
background, so that a background originating from the ra- 
dio photosphere would, on average, lead to a brighter set of 
masers. 

An interesting possibility is that it is the background 
which breaks the symmetry of the maser shells in real stars, 
leading to the distribution of spots that we see, instead of 
an unbroken ring. In the old model, we appealed to vari- 
ations in the SiO abundance to break the symmetry, but 
as we now have this as a smooth function of radius (see 
Fig. [3]) it is less acceptable to do this. A purely specula- 
tive idea at present is that the maser sites lie above regions 
which are like extensions to the radio photosphere, with en- 
hanced electron abundance from above average ionization of 
the alkali metals, and consequently high input backgrounds 
to overlying masers. These sites would then have the advan- 
tage over surrounding regions of needing shorter amplifying 
columns to reach the same maser intensity. Asymmetries 
of this type in the molecular layer, which forms the op- 
tically thick band radius at 8 jj^m have a lready been found 
("Weigel t et aIlll996l : [Ragland et al.ll2008h . although it is not 
clear that these would extend, in wavelength, to the radio re- 
gion. The choice of maser sites in our comp utations remains 
one of random selection (^Humphrevs et allflQQd ). 

3.3 Pointlike Masers 

The use of the LVG approximation for the transport of the 
pumping radiation simplifies the radiation transfer problem, 
but there is a considerable price to pay for this simplicity. 
The approximation is essentially local, so that a region con- 
sitituting the maser is effectively defined by the distance 
(which can vary with direction) required for the radiation 
Uneshape to be shifted out of resonance with the response 
lineshape of the molecules. When we pick the coordinates 
of a maser site, there is still an assumption that the physi- 
cal conditions within this, velocity bounded, LVG zone are 
equal to those at the selected coordinates. 

3.4 LTE Chemistry 

The abundances of the atomic and molecular species used 
in the current model are based on LTE chemistry. This 
amounts to an assumption that the network of chemical 
reactions comes to equilibrium on a timescale which is 
short compared with important hydrodynamic timescales 
in the envelope. It is likely that non-LTE chemistry is im- 
portant is establishing the time-varying populations of a 



number of species, inc luding CO, H CN,C S and SiO, as pre- 
dicted in models by ICherchnefl (|2006l V Detection of all 
these molecules in a sample of AGB stars of widely vary- 
ing C/0 abundanc e strongly supports the model predictions 
(jPecin et aLlbOOSh . There is a brief discussion of non-LTE 
effects pertaining to m odels of the type used in the present 
work in [Ireland etaP (jioQS ). 



4 COMPUTATIONS 

4.1 IR Pumping Radiation 

We briefly consider the conversion of the dust optical effi- 
ciencies described in Section 12.21 to the inputs required by 
the maser model: namely effective grain areas and angle- 
averaged intensities (J) at the wavelength of every electric- 
dipole-allowed transition in the SiO model. The areas were 
trivially obtained from the optical efficiencies, since the 
grains are modelled as having a unique radius. 

To compute the angle- averaged intensities, we solve, 
for each tabulated dust wavelength and radius in the dust 
model, the radiative transfer equation, 

dr 1 + (ks/ka) + {ct/ka) 

l + {KA/Ks) + {a/Ks) l + K/a' ^' 

which is written for the most complicated situation consid- 
ered, with a two-component dust shell of AI2O3 (subscript 
A) and silicate (subscript S). Absorption coefficients (k,) 
and scattering coefficients (a) are written without a sub- 
script when summing over both dust types, for example 
cr = (JA-\-(Js- The optical depth is given by the standard ex- 
pression, dr = {Hi-\-a)ds, where ds is the element of distance 
along a ray. With the single-radius dust approximation, it 
is possible to write the ratios that appear in eq.{[T]) directly 
in terms of the optical efficiencies for absorption, Qabs, and 
scattering, Qsca, for example 

K,S _ nsQabs,S /^N 
KA TlA^ahs^A 

where ns {tia) is the number density of silicate (AI2O3) 
dust. Similar expressions can be used for the other ratios. 
The thermal emission terms contain Planck functions, B{T), 
for the appropriate dust temperature, and the scattering 
term contains the angle- averaged intensity, J, defined in the 
spherical geometry used here as 

We note that as eq.([T|) contains J, the solution is iterative, 
with specific intensities for each value of /x = cos^ being 
computed with the most recent estimate of J, before a new 
value is found from eq.©. The boundary conditions used in 
the solution of eq.([T]) were to use a Planck function at the 
kinetic temperature of the appropriate SiO band radius for 
values of /i that intersected such an optically thick layer. For 
other values of /x, a Planck function at the CMB temperature 
was used. A value of J was calculated for each radius in the 
hydrodynamic solution, and for each dust wavelength. 

We next calculate the dust contribution to the pump- 
ing rates in each SiO transition. For a transition with upper 
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Figure 6. A maser site shown in the cavity between an optically 
thick layer, marked as the star, and the surrounding dust shell 
(case (ii) in the text). Rays at various possible angles are shown 
converging on the maser site, where they are angle averaged. The 
size of the optically thick layer is different, in general, for each 
SiO pumping band and pulsational phase. 



level, and lower level, /, the dust contributions are Bui J 
and BiuJ, where Bui and Biu are the Einstein B- coefficients 
for the transition. We first assign each transition to either 
one of the IR pumping bands, or as 'radio', a blanket des- 
ignation for all the pure rotational lines. The input data 
for the dust routine is then read in: the dust areas tab- 
ulated by wavelength, the dust temperatures and number 
densities tabulated by radius, and the dust mean intensi- 
ties, computed as discussed above, tabulated by both radius 
and wavelength. These mean intensities are required for cal- 
culating the scattering contribution to the dust radiation 
field. Then, considering the radius of a single maser site, a 
geometrical dilution factor is calculated relative to the op- 
tically thick band radius at each pumping band. In a loop 
over the SiO spectral lines, we then classify the maser site 
as in one of three basic positions for the current line. Case 

(i) is that the site is inside the band radius, in which case 
the pump in the current line corresponds to a black-body 
at the temperature of the site radius and wavelength of the 
line; there is no dust contribution. Case (ii) is that the maser 
site is in a dust-free gap between the band radius and the 
dust shell, for example Fig.O The final possibility, case (iii), 
is that the maser site is inside the dust shell. In both case 

(ii) and case (iii), a radiative transfer equation is solved for 
the current transition. Rays were traced at a number of an- 
gles through the envelope, with each ray terminating on the 
maser site. Values of the dust areas, and mean intensities 
at the line wavelength were obtained from spline interpola- 
tions from the input tables for each site radius. A given ray 
had two possible starting conditions: an optically thick layer, 
yielding a black-body at the temperature of that layer, or an 
interstellar value, taken to be the CMB. Case (iii) required 
the stellar contribution to undergo some extinction due to 
the dust. Angle- averaging of the rays produced the required 
mean intensity at the wavelength of a particular line, and 
this was trivially converted to a pump rate. 



Table 3. Jobs run to investigate the effect of variation of the dust 
shell only. All other parameters are for the model with structure 
M21n (ilreland et al.ll2004br) . 



Job Number 


Structure 


Dust Epoch 


1 


M21n 


D 


2 


M21n 


C 


3 


M21n 


B 


4 


M21n 


A 



4.2 Size of Radio Photosphere 

We have noted in Section 13.11 that there is a considerable 
uncertainty in the size of the radio photosphere for the maser 
lines. We here show the difference between a model based on 
the standard value, which is equal to the 8 /xm radius, and 
an otherwise identical model in which the radio photosphere 
has been taken to coincide with the IR photosphere. In both 
cases we use Model 1 (see TablefS}, and the same distribution 
of maser spots. In the first version of the model (la) we allow 
masers with nominal positions outside an IR photospheric 
radius of 1.21 Rp (the 1.04 /im value) to grow; in Model lb, 
we exclude masers inside 2.20 Rp, corresponding to the 8 /am 
radius. This exclusion is based on the idea that, though SiO 
molecules may exist at the site, and these may be sucessfully 
pumped, optically thick layers at larger radii, and at the 
maser frequency, will prevent external observers from seeing 
the emission. 



4.3 Variation of dust envelopes 

Initially, we considered a single hydrodynamic solution, 
M21n, and ran the maser code for 1500 maser sites, as in the 
old model, for each dust epoch in turn. We note that M21n 
is actually paired with dust epoch D, but we ran it with the 
other three epochs simply to look for the effects of varying 
the dust component of the radiative pump. The models test- 
ing the variation of the dust envelope are listed in Table (3] 
The seed for the random number generator was the same 
for all these models. As they also have the same structure, 
which means the same parent radius for the star, the maser 
positions were the same for all the models. Results for the 
dust variation models appear in Section [521 



4.4 Phase Samples 

In this set of computations we ran jobs for four hydrody- 
namic solutions, each with its standard dust epoch. Details 
of these jobs appear in Table |4l Note that job 1 appears in 
both the dust variation and phase sample data. The mod- 
els all used 1500 maser sites, randomly distributed between 
the (phase dependent) IR photospheric radius (see Table 2]) 
and an outer radius of 5Rp. In view of the small number 
of phase samples, we did not atter apt to advance the mode l 
with phase, the method adopted in'Humph revs et al.l (|2002[ ), 
but instead we re-randomized the maser positions for each 
new model. Results appear in Section [531 
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Table 4. Jobs run to investigate the effect of stellar phase on the 
intensity an d distribut ion of maser spots. Structure designations 
are from llreland et al.l (2004b). Phases are based on the optical 
light curve. Stellar radii (Column 5) are based on the IR photo- 
sphere, at 1.04 ia.m, and given in units of the 'parent' radius of the 
star, where Rp = 1.809 x 10-*^^ cm. 8.13 fxmradii, assumed equal 
to the size of the radio photosphere, are also in parent units. 



Job 


Structure 


Dust 


Phase 


Ri.04/Rp 


Rs.i3/Rp 


1 


M21n 


D 


0.10 


1.21 


2.20 


5 


M22 


A 


0.25 


1.10 


1.77 


6 


M23n 


C 


0.30 


1.03 


1.73 


7 


M24n 


B 


0.40 


0.87 


1.69 




Velocity (km/s) 

Figure 7. Spectra of the three SiO maser lines ^' = 1,J = 1 — 
(solid line), 'u = l,J = 2 — 1 (dashed line) and 'U = 2,J = 1 — 
(dotted line) for Model la (top) and Model lb (bottom). These 
model differ only in the assumed size of the radio photosphere, 
which is smaller in Model la (upper graphs). 

5 RESULTS 

5.1 Size of radio photosphere 

In Fig. [71 we show spectra from Model la (top) and Model lb 
(bottom). Recall that in Model la, the radio photosphere is 
small, coinciding with the IR photosphere, which at this 
phase has a radius of 1.21 Rp. By contrast, in Model lb, the 
radio photosphere coincides with the 8.13 iJ,m radius, equal 
to 2.20 Rp. The three spectra shown in each panel are the 
three most commonly observed lines: v — l^J— 1 — and 
J = 2 — 1 and ^7 = 2,7 = 1 — 0. In both lines from v — 1^ 
Model la (top panel) produces extremely high amplification 
masers. These are emitted by a population of maser sites, 
just outside the small radio photosphere of this model. We 
refer to this population as the 'bright, inner object' to distin- 
guish them from more modest maser emitting sites at larger 
radii in the envelope. 

The V — 2 transition has the same shape in both panels, 
so the bright, inner, objects do not emit significantly in this 
line. Emission from the bright inner masers in both the v — 1 
transitions can be seen to be blue-shifted by ^5kms~^ rel- 
ative to the emission from the outer population. This shift 
arises naturally if the inner population of masers is associ- 
ated with post-shock gas inside the inner shock (see Fig. [1]) 

In Fig. [8] we show synthetic maps of the 43 GHz v — 




-4-3-2-101234 
R (parent units) 

Figure 8. Modelled maps of the v = l,J=l — transition with 
radio photospheric radii of 1.21 Rp (top) and 2.20 Rp (bottom). 
In the top panel, the log to the base 10 of the maser amplification 
is clipped at a value of 3, for any individual spot, to render the 
outer ring visible. 

1, J = 1 — masers, the most commonly mapped transition 
with VLBI instruments. In the upper figure, we again have 
the smaller radio photosphere, allowing the appearance of 
an inner ring of very bright masers, with a projected radius 
close to ^IRp. The lower panel shows the maser distribu- 
tion with the larger radio photospheric radius of 2.20 Rp. 
The maser ring in the lower panel has properties which are 
similar to those of a 'quiet' phase in the old model, and a 
radius which is in good agreement with observations. If the 
spots of the inner ring had similar dimensions and beam an- 
gles to those at larger radii, they would be a factor of 10^ 
times brighter than those found in the old model, or at larger 
radii in the current model. As masers of this power are not 
observed, we conclude that such a small radio photosphere is 
not reaUstic, and we adopt the larger, standard, value from 
here on. We note that, even if we adopt the smaller radius for 
the radio photosphere, the masers still form in two distinct 
rings: there is not a population of maser spots of similar 
brightness extending from the outer ring inwards towards 
the inner group, though the mapping software used to gen- 
erate Fig. [8] suggests that there may be. In fact, if we adopt 
a radio photospheric radius of 1.39 Rp, corresponding to the 
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2.03 ^im pumping band, the spectra and maps are almost 
identical to Model lb, which features only the outer ring. 
This intermediate radius is already large enough to exclude 
all the very bright spots with no observational counterparts. 

An interesting feature of the outer ring (Fig. [8]) is 
a 'missing' sector in the (observational) South- West. This 
arises purely from random selection of sites, but such gaps 
do sometimes appear in the maser rings of real stars. 

We adopt the larger radio photosphere from now on, 
corresponding to the S.lS^m pumping band. Therefore, the 
extremely bright, inner, ring of masers is excluded. How- 
ever, the fact that this inner ring can be produced raises an 
interesting point about the origin of the maser sites. If we 
assume that the motion of a typical maser spot is outward, 
more than inward, then we expect new spots to feed into the 
distribution near its inner edge, and therefore close to the 
radio photosphere. How do these new spots appear, and do 
they predominantly 'switch on' at a particular phase? The 
old model did not address this question because it set an ini- 
tial random distribution of spots at model phase zero (when 
a new shock left the stellar photosphere) and subsequently 
adjusted the positions of these spots with phase, following 
the velocity field in the envelope. There was no source of 
new spots. The new model may be able to derive a phase 
for the appearance of a new ring once more phase samples 
are available. 

5.2 Variation of dust envelopes 

We find that the effects of varying the dust epoch for the 
single hydrodynamic solution, M21n, vary from subtle shifts 
to profound change, depending on the transition in question. 
As an example, we show in Fig. (9] model spectra for three 
commonly observed lines. Beginning with the top panel, the 
43 GHz V — 1,J = 1—0 line, we see that the spectra have 
very similar shapes for all four epochs, with a fairly subtle 
change in the amount of amplification introduced by vary- 
ing the dust radiation field. In this line, Model 3, with dust 
epoch B, produces the brightest spectrum, with a peak am- 
plification factor of 108, followed by Model 1 (epoch D dust) 
with 84, Model 2 (epoch C dust) with 80, and finally Model 4 
(epoch A dust) with 41. Varying the dust radiation field 
therefore changes the peak maser amplification by a maxi- 
mum factor of just over 2 in this line. The v — 2^J— 1 — O 
line (bottom panel in Fig|9]) also responds quite subtly to 
changes in the dust envelope, with very similar shapes for 
all four epochs. However the brightness order is different 
from the analogous line in v — 1: here Model 1 yields the 
strongest line, though Model 4 is still the weakest. 

In the 86 GHz line (t' = 1, J = 2 — 1) we see a more 
profound, qualitative, effect of varying the parameters of 
the dust envelope. Spectra for this transition appear in the 
middle panel of Fig. (9] Model 4, with dust epoch A, which 
yields the weakest spectrum in both 43 GHz lines, introduces 
a strong additional peak on the blue side of the spectrum at 
86 GHz, which is not present at the other 3 epochs. 

A more general view of the effect of the dust radiation 
field on the SiO masers is presented in Fig jlOl Here we plot 
the peak amplification factor found in each spectrum. Each 
panel of the figure is a different vibrational state. Different 
symbols denote the dust epoch, and peak values are marked 
against upper J-quantum number for transitions J = 1 — 
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Figure 9. Model spectra for the maser transitions 'u = l,J=l— 
(top panel), 'u = l,J = 2 — 1 (middle) and 'U = 2,J = 1— (bottom 
panel). Each panel has four versions of the spectrum correspond- 
ing to: Model 1, dust epoch D, (solid line), Model 2, dust epoch C, 
(long dashes), Model 3, dust epoch B, (short dashes) and Model 4, 
dust epoch A (dotted line). The hydrodynamic solution is M21n 
in all cases. 



up to J = 10 — 9. The main point to be drawn from Fig. [10] 
is that there is something special about Model 4, which uses 
dust epoch A, and is marked by the symbol □. Model 4 is 
the only one of the four to generate any masers in = 0, 
where inversions are only generated in transitions where the 
upper rotational level has even J. In the other vibrational 
states, only Model 4 has masers for upper J > 5; at lower J, 
the peak amplification values generated by the other three 
models are closely grouped, while the value for Model 4 is 
often very different. The middle panel of Fig|9]is an example 
of this. The even-odd pattern produced by Model 4 in = 
is also echoed through the other ^;-states: masers with upper 
J even are favoured in v = 0, 1, but masers with upper J odd 
are favoured \n v = 2 and 3. We cannot, at present, explain 
the pecuharity of Model 4 in detail, but we note that in 
FigEl the Epoch A dust has the second-smallest inner dust 
shell radius, but a significantly higher dust temperature at 
its inner radius than the other three models. 

If the inner ring of masers is considered, the effect of 
changing the dust parameters was much more marginal. 
The weaker, but 'standard' outer ring of spots are much 
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Figure 10. Peak amplification factors found in the spectra of 
the first ten rotational transitions (upper J-quantum number on 
the X-axis) of the first four vibrational states of SiO (panels as 
marked with v increasing from bottom to top. Values for different 
dust epochs are marked with a different symbol: + for Model 1 
(Epoch D), X for Model 2 (Epoch C), * for Model 3 (Epoch B) 
and □ for Model 4 (Epoch A). 
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Figure 11. Spectra of the masers found in a dust-free model, 
based on the hydrodynamic solution M21n. Rotational transitions 
have the following line styles: J = 1 — (solid), J = 2 — 1 (dashed), 
J = 3 - 2 (dotted), J = 4 - 3 (chained). 



line was excluded becasue of its interaction with a water 
transition. 

When we compare the dust-free models to those with 
dust shells, we find that the largest effect is in ^; = 1, = 
1 — 0. The shape of the spectrum is similar, but if we compare 
the solid line in the top panel of Fig. [11] to the lines in the 
top panel of Fig. [9l we see that the peak in the dust-free 
case is brighter, by a factor of approximately 2, than the 
brightest case with a dust shell. For the other lines in Fig. [9] 
[v — 2^ J — 1 — ^ and ^7 = 1,7 = 2 — 1) we again find similar 
shapes (with the exception of dust epoch A) and smaller 
changes in maser amplification. Overall then, the presence 
of AI2O3 dust can have both positive and (mildly) negative 
effects on maser pumping. We intend to investigate the effect 
of silicate dust, if present, in future work. 



more strongly affected individually, and when grouped into 
a spectrum, than those which were very intense. 

To test the 'absolute' effect of dust, an additional model 
was run without a dust shell. In this model, pumping is 
entirely due to the radiation from the star (in Model M21n) 
and kinetic collisions. We show spectra for those lines in 
V — 1^ V — 2 and v — 3 which showed significant maser 
emission in Fig. 1111 As elsewhere, the v = 2, J = 2 — 1 



5.3 Phase Samples 

In the four phase samples listed in Table HJ no masers were 
found in V = 0. In V — 3, weak masers were found in rota- 
tional states up to J = 5 — 4, with the greatest amplification 
factor, of 19.2, at phase 0.1 (Model 1) in v = 3, J = 2 - 1. 
All the V = 3 masers followed a pattern of producing the 
strongest maser emission at phase 0.1, with a monotonic de- 
cay to phase 0.4. Only J = 1 — and 2 — 1 masers survived 




Figure 12. Spectra of the bright masers in = 1 for the phase 
samples in Table 0] Each panel is for the transiton as marked. 
Spectra in each panel have the following line styles: (f)=0.1, 
Model 1, (solid hne), (j)=0.25, Model 5, (long dashes), (f)=0.30, 
Model 6, (short dashes) and (/)=0.40, Model 7, (dotted hne). 

to reach phase 0.4. In Fig. [12] we plot spectra of the three 
masers from v — 1 which achieve amplification factors of 
greater than 20. These form in the lowest three rotational 
transitions. Figure [13] shows the spectra of bright masers 
from V — 2. Here, we ignore the J = 2 — 1 line because 
of its well-known interaction with a water transition, so the 
plot consists of the masers in the J=l— 0, J — 3 — 2 and 
J = 4 — 3 transitions. 

The main points to note from Fig. [12] are (i) the J = 
1 — line is the most stable: it is the only line which still 
emits as a maser at phase 0.4; (ii) the most favourable phase 
for masers is 0.25 (Model 8): both the J = 1-0 and J = 2-1 
lines are strongest at this phase, and in J = 3 — 2 the maser 
from the next phase is only marginally stronger; (iii) in the 
J = 1 — line, there appears to be a systematic redward 
shift of the spectral peak with phase, which does not seem 
to be the case in the other two lines. 

In Fig. [131 we find that, as for ^; = 1, the v = 2, J = 
1-0 line produces the most stable spectrum, with the only 
surviving emission at phase 0.4. The redward shift of the 
spectral peak with phase is apparent for all three lines shown 
in Fig. I13[ with a very distinct shift between phases 0.1 and 
0.25 in the J = 1 — and J = 3 — 2 transitions. This shift 
significantly changes the shape of the spectra, while the peak 



Figure 13. As for Fig. [12] but for v = 2. 

moves to the red by ^7kms~^. The brightest phase is 0.25 
for the J = 1 — line, but 0.1 for the other two. 

We now consider the spatial distribution of maser spots 
as a function of stellar phase in four transitions: J — 1 — in 

V = 1 and ^; = 2, V = 1, J = 2 - 1 and ^; = 2, J = 3 - 2. The 
v = l,J=l— maser shell is shown in Fig. 1141 The main 
change between phase 0.1 and the brightest phase (0.25) 
is a decrease in the ring radius, as judged by the bright- 
est spots. After phase 0.25, there is a fading of the masers, 
with only a single obvious object remaining in the maps for 
phase 0.4. For the other 43 GHz line, v = 2, J = 1 — 0, we 
show its evolution with phase in Fig. 1151 The behaviour is 
broadly similar to that of its counterpart in v = 1, but it is 
clearer that, in the v — 2 transition, something catastrophic 
happens to the maser distribution between phases 0.1 and 
0.25. Although the masers are brighter in the latter phase, 
the ring is smaller, indistinct, and based on just a few very 
bright objects. These fade almost to extinction by phase 0.4. 
It is noticeable that in the map at phase 0.1, the region cov- 
ered by weak maser emission is larger than for v — 1, so the 

V = 2 maps are plotted out to 5, rather than 4Rp. The ap- 
parent North- South line of symmetry in the later phases in 
Fig. [15] results only from the random spot distribution and 
the mapping software, and does not result from any feature 
of the hydrodynamic solutions or dust shell. A comparison 
of Fig. [14] and Fig. [15] phase by phase, suggests that, for 
the bright masers at least, the v = 2 spots occupy a smaller 
ring than those in v = 1. In Fig. 1161 we show the maser 
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_4 -3-2-101254 

R [parent onits) Figure 15. As Fig. [H but for the V = 2, J = 1 - transition. 



Figure 14. Maps of the spot distribution for masers in the v = 
1, J = 1 — transition for the four phase samples, as marked. The 
colour scale is in the logarithm to the base 10 of the amplification. 



shell of the 86 GHz transition, ^7 = 1,7 = 2 — 1. Only three 
phases are shown, as there is no maser emission at phase 
0.4. In this line, the ring structure is distinct through the 
three phases shown, but, as with the 43 GHz lines, there 
is evidence of a distinct change in structure, with a mod- 
estly strong, large radius maser ring (~3Rp) replaced by 
an intense ring of radius <2Rp. The final set of maps we 
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Figure 16. As Fig. [M] but for the v = l,J = 2 — 1 transition. 
No map is shown for phase 0.4 because the model produced no 
emission in this hne for this phase. 



Figure 17. As Fig. [TH but for the = 2, J = 3 - 2 transition. 
No map is shown for phase 0.4 because the model produced no 
emission in this line for this phase. 



show, which is not yet observable with a VLBI interferom- 
eter, is the V = 2, J = 3 — 2 transition in Fig. 1171 Like 
V — 1,J = 2 — 1, there is no emission at phase 0.4. The 
first two phases have a very similar peak amplification, as 
seen in the spectrum in Fig. I13[ but the spatial distribution 
of spots has still suffered a similar change in distribution to 
the other mapped transitions between phases 0.1 and 0.25. 
This change is also reflected in the spectrum of this line. 
The maser ring is arguably less rich than those from lower 
rotational states. 



Extended structure visible in the all the maps in this 
work (Figs. ll4]fT71 and Fig. [8]) should be viewed with great 
caution. All the masers are actually pointlike (see Sec- 
tion [33]). The GNUPLOT software used to make the maps 
includes a gridding algorithm which smooths the spot distri- 
bution into pixel values on a regular grid. The colour palette 
is then applied to these pixel values. However, there are typ- 
ically several hundred low amplification maser spots which, 
when gridded and smoothed by the mapping software, yield 
the extended structure. If our distribution of low amplifi- 
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cation points generate an equivalent flux to weak extended 
emission in real stars, could they explain observed discrepan- 
cies between the integrated flux from VLBI images and the 
flux from autocorrelation spectra? If we simply divide our 
spots into two classes of bright and weak, and assume that 
the masers in both classes have the same sizes and beam- 
ing angles, then the ratio of fluxes produced by the bright 
and weak objects is NbIb/{NwIw), where N is the number 
of objects, subscripted with b for bright, or w for weak, and 
similarly for brightnesses, or specific intensities, /. As an 
example, we take the map for phase 0.1 in Fig. [M] and set 
the dividing line at spots with amplification factors of 200. 
We then have 6 bright and 256 weak spots. To obtain sim- 
ilar fluxes from the weak and bright populations we then 
need an average weak spot to be no fainter than 1/43 of 
the brightness of an average bright spot. On this basis our 
distribution of weak spots might possibly provide the flux 
missing from VLBI images. 



6 DISCUSSION 

The radio photosphere plays a vital role in reconciling the 
new model with observations. By selecting a rather large 
value of this parameter, a radius similar to 2Rp, a group 
of extremely bright masers in an inner ring of projected ra- 
dius ~lRp is excluded, leaving the remaining masers in a 
ring which is consistent with both observations and the old 
model. It is also useful to look at what else is obscured by se- 
lecting the large radio photosphere. In particular, the strong 
inner shock with a velocity change of >20kms~^, which is 
shown in Fig.[T]at a radius of 2.9 x 10^^ cm, would not be vis- 
ible in the radio. Our radio photosphere at this phase (0.1) 
has a radius of 2.2 Rp, or 3.98 x 10^^ cm, so we would only be 
able to see the outer, highly damped, shock in Fig.[T] and the 
maximum observable velocity range in the envelope at radio 
frequencies is only ^7kms~^. Velocities of this magnitude, 
but not the 20km s~^ of the inner sho ck, are in agreement 
with observations (Reid & Menten 2007), which measure the 
maximum velocity from the proper motions of SiO masers, 
and link these spatially to the radio photosphere with VLA 
continuum measurements. The damping of the shocks in the 
hydrodynamic solutions used in the current work is there- 
fore an important improvement over the old model, in which 
the shocks were less well damped, and propagated at higher 
speeds for a given radius in the envelope. 

The effects of dust can be both constructive and detri- 
mental to the maser pump. The most spectacular effect is 
that the dust in Model 4 (epoch A) has a qualitatively dif- 
ferent effect on maser emission when compared to the other 
three dust epochs, used in conjunction with the same hydro- 
dynamic solution (M21n). We investigate briefly the origin 
of this behaviour, for which the v — 1, J = 2 — 1 line is a good 
example (see Fig. [9]). On close examination, a single bright 
model spot is responsible for much of the extra peak in the 
spectrum of this line. This spot is positioned at a radius 
(unprojected) of 3.59 Rp, which is rather large, but the ring 
radius for this line is also larger than average (see discus- 
sion below and Fig. [16]). This radius places the bright spot 
well within the epoch A dust shell, very close to the 9th 
shell radius. We compared the dust mean intensity at the 
spot radius for all four epochs, and all the SiO IR pumping 



bands, but there is nothing peculiar about epoch A in this 
respect. It is therefore not the direct pumping power of the 
dust that makes epoch A different. A feature of the epoch A 
shell which is different from the other three epochs, with re- 
gard to the selected spot position, is a combination of high 
dust temperature, density and considerable depth into the 
shell. Epoch D has a higher dust temperature (1070 K versus 
930 K), but the spot lies only between the 3rd and 4th radii. 
The spot lies deeper inside the dust (near the 11th radius) 
in epoch B, but this shell has lower density and tempera- 
ture (715 K). The most likely explanation for the particular 
behaviour of masers with dust from epoch A is therefore an 
optical depth effect, with the dust modifying the stellar IR 
continuum in the pumping bands. We intend to investigate 
the effects of dust further in future work. 

The phase information in the present work is rather 
limited, so we have not attempted to follow the evolution of 
maser rings through the stellar cycle in a consistent way. We 
note also that in the modelling of S Ori, which used the the 
same hydrodynamic solut ions and dust shells, the fits to the 
observed optical phases (|Wittkowski et al ]|2007'), and see 
Table [2l have a sequence which does not match that of the 
phases drawn directly from the theory. In spite of this, and 
following the theory-based phases, we find that the maser 
peak for all lines is in the phase range 0.1 to 0.25. This 
is consistent with the maser peak lagging the optical light 
somewhat. Given the error of 0.1 in the absolute phases from 
the model, this phase lag is consistent with typical observed 
lags between SiO masers and the optical light curve. 

A distinct feature of all the maps is a structural change 
between phases 0.1 and 0.25, where the maser ring shows a 
decrease in radius and, usually, an increase in spectral out- 
put. This change in the structure of the map is accompanied, 
in some lines, by a redward shift in the spectral peak. There 
is also a shift from a large number of maser spots of modest 
brightness, to a smaller number of objects, a few of which 
are very intense. This behaviour may have some underlying 
cause related to the envelope dynamics, but as spot posi- 
tions have been re-randomised, and not evolved, between 
the two phases, we can only suggest that this is an inter- 
esting feature, worthy of further investigation, rather than 
a hard result. 

The new model allows us to be considerably more pre- 
cise than the old model in studying the spatial relationship 
of the maser rings to each other, and to other important 
structures in the model, such as shocks, dust shells and op- 
tically thick layers. We display this information in Table [S] 
The first six columns of this table display model parameters, 
which are independent of the maser model. The remaining 
6 columns show maser ring radii for the four transitions 
mapped in Fig. 11411171 Of these six columns, the first four 
are ring radii from the present model, whilst the final two are 
from observations of S Ori, matched to theoretical phases as 
in IWittkowski et al.1 (|2007). All radii given in Table [5] are 
relative to the IR photosphere (at 1.04 iJ,m) for the particu - 
lar phase for easy comparison with lWittkowski et al.l (|2007f ) , 
Table 7 and Fig. 12. 

Without even considering any maser ring radii, we can 
see that the 8.13^Lm band radius (Column 6) is closely re- 
lated to the position of the innermost shock (Column 3). 
The radii are the same to 0.1Ri.o4 fo three epochs, but 
markedly different for epoch D (phase 0.1). This finding 
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Table 5. Key radii for maser shells, and their association with shocks and optically thick band radii. All 
radii are written in terms of the IR photospheric value, i^i.o4, at the appropriate phas e. Columns are as 
follow s: 1 - the hydrodynamic solution and dust epoch (the same as the observing epoch in I Wittkowski et al.l 
(20071)): 2 - the stellar phase from theory (optical scale); 3 - radius of the inner shock; 4 - radius of the outer 
shock (not apparent at phase 0.4); 5 - inner radius of the dust shell; 6 - optically thick band radius in the 
Av = 1 SiO pumping band; 7 to 10 - maser ring radii (in sky projection) with sample stan dard deviations for 
the tr ansitions marked from the present model; 11 and 12 - observed maser ring radii from I Wittkowski et al.l 
(|2007h for transitions as marked. In columns 3 and 4, we have listed the radii for the same shocks. At phase 
0.4, the original outer shock has damped out, so we have left the entry for Rs2 blank at this phase. A new 
shock has formed at phase 0.4, with a radius of I.OR1.04, inside the original inner shock, listed as Rsi- 



Model Parameters Maser Ring Radii (photospheric units) 



Model+dust 






Rs2 




Rs.is 


R43.1 


R42.8 


^v^l,J=2-l 


R7; = 2, J = 3-2 


■p 06s 
-^43.1 


■p obs 
-^42.8 


M21n+D 


0.1 


1.3 


2.5 


2.4 


1.8 


1.9±0.4 


1.9±0.4 


2.4±0.7 


2.0±0.4 






M22+A 


0.25 


1.6 


2.6 


1.8 


1.6 


1.9±0.3 


1.8±0.3 


1.8±0.5 


1.8±0.3 


2.2±0.3 


2.1±0.2 


M23n+C 


0.3 


1.7 


2.8 


2.2 


1.7 


1.9±0.2 


1.8±0.2 


1.9±0.2 


1.9±0.2 


2.1±0.3 


1.9±0.2 


M24n+B 


0.4 


2.0 




2.0 


2.0 


2.1±0.5 


2.4±0.4 






2.4±0.3 


2.3±0.4 



strongly suggests that the whole way of thinking about 
shocks propagating vast distances from some almost-fixed 
stellar photosphere, which caused many problems for the 
old model is simply outdated. For a significant portion of 
the cycle at least, the S.lS^Lmradius is tied to a shock, so 
that the delay between irradiation and shock impact would 
be 25 d, or less. Consequently, many arguments which elim- 
inate some element of collisional pumping on the grounds of 
long shock propagation times, relative to a distant, detached, 
stellar photosphere, are not relevant to resolving the pump- 
ing debate. The physical changes in the maser shells between 
phases 0.1 and 0.25 may be caused by the shift from a state 
where the 8.13 ^tm radius and inner shock are detached, to 
one where they are closely coupled. 

The outer shock appears to be well outside the maser 
zone, except at phase 0.1, where it may be associated with 
the large ^7 = 1,7 = 2 — 1 ring (see Column 9). At most 
phases, the inner shock, the 8.13 ^tm radius and the maser 
ring radii of all the lines are within 0.3Ri.o4 of each other. 
This association does not involve the inner shock at phase 
0.1. The inner edge of the AI2O3 dust, Rin (Column 5) is less 
well associated with masers, inner shock and 8.13 iJ,m radius. 
It is close to the maser ring radius at phase 0.25 and 0.4, 
but not at phases 0.1 and 0.3. 

The maser rings themselves are closely associated for 
almost all transitions and phases. The exception is the large 
ring at phase 0.1 for the v = l,J = 2— 1 transition. The ring 
radius for this transition is not peculiar at the other phases. 
In the first three phases, the mean ring radius of the v = 
2, J = 1 — masers is the same, or slightly smaller than that 
for the analogous transition in v = 1. This trend is swapped 
for phase 0.4, but this phase has a large standard deviation, 
due to the small number of surviving spots, and a rather 
poorly defined ring. At this phase, close to optical minimum, 
the masers from both J = 1 — (43 GHz) transitions have 
their largest mean ring radii in both observations and theory. 
The model is therefore consistent with previous models and 
with VLBI observations which show that the v — 2^ J — 
1 — ring is usually slightly smaller. The large ring for the 
V — 1^ J — 2 — 1 masers at phase 0.1 is not in agreement with 
the old model, but is m ore in line with VLBI observations 
(jSoria-Ruiz et al.ll2007l ). The modelled ring radii at 43 GHz 
are also smaller than their observational counterparts, for 



the three phases where observational data exist. However, 
the ring sizes are consistent, given the standard deviations 
attached to these ring radii. 

Maser spectra in the 43 GHz (J = 1 — 0) lines have 
velocity extents of ~10kms~^, a value which is consistent 
with the observational s pectra of these lines towards S Ori 
(jWittkowski et al.ll2007l ). The LSR velocity of S Ori is not 
very well known, so we have not computed Doppler biases 
for our spectra, though these model spectra are certainly 
concentrated to the red side of the velocity range. 

Perhaps the largest drawback of the current work, in 
terms of matching VLB A observations of S Ori, is in the vari- 
ation of spectral intensity. In the three observed phases, the 
ratio of the spectral maxima for the brightest and dimmest 
phases is about 1.5 for the ^7 = 1,7=1- line, and about 
2 for 'i; = 2, J = 1 — 0. By contrast, the fading of the model 
masers towards phase 0.4 gives analogous ratios of approxi- 
mately 70 and 1000. However, as has been noted previously, 
our models are not evolved in phase, but use repositioned 
spots, so the numbers could be very different if the masers 
followed the gas motions. There are also only a small num- 
ber of bright spots, so these ratios are very vulnerable to 
statistical fluctuations. 



7 CONCLUSIONS 

The new combined maser and hydrodynamic solutions pro- 
duce maser emission in rings of typical size 2.2 IR photo- 
spheric radii. The ring structures are similar in many re- 
spects to those produced by the old model. An analysis of 
ring radii confirms the usual observational result, and the 
prediction of the old model, that the v — 2 ring is smaller 
than the v — 1 ring at 43 GHz, whilst a large ring is pos- 
sible at some phases in the v — l^J — 2 — 1 transition at 
86 GHz. The overall range of maser ring radii, from 1.8-2.4 
IR photospheric radii, for the 43 GHz transitions in v — 1 
and V — 2 is also consistent with VLBI experiments where 
maser rings and photospheric radii have been compared at 
the same phase (|Boboltz k, Wittkowskj l2QQ5l : iFedele et al.1 
I2OO5I : IWittkowski et al.ll2007l ). 

If a radio photosphere of approximately 2 IR photo- 
spheric radii is used, observable shock velocities are in agree- 
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ment with actual radio continuum observations, as high ve- 
locity shocks, such as the inner shock at phase 0.1, are hid- 
den. A radio photosphere of this size also prevents the ap- 
pearance of extremely bright masers very close to the IR 
photosphere. 

Dust can both suppress and enhance maser emission, 
and there is definitely a radiative component to the pump- 
ing scheme. However, the close association of shock and 
8.13 micron radii suggests that collisional and radiative 
pumping are closely associated spatially, and therefore tem- 
porally. 

Further work is required to consider the effects of dust 
in detail, and to produce a fully phase-sampled model of one 
or more complete stellar cycles. 
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